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Abstract: We detail the continued fraction expansion of the square root of 
a monic polynomials of even degree. We note that each step of the expansion 
corresponds to addition of the divisor at infinity, and interpret the data yielded 
by the general expansion. In the quartic and sextic cases we observe explicitly 
that the parameters appearing in the continued fraction expansion yield integer 
sequences defined by bilinear relations instancing sequences of Somos type. 

The sequence 3, 2, 1, 1, 1, 1, 1, 2, 3, 5, 11, 37, 83, ... is produced by 
the recursive definition 

Bh+3 = (Bh-iBh+2 + BhB] l+ i) I Bh-2 (1) 

and consists entirely of integers. 

On this matter Don Zagier comments |26| that 'the proof comes from the theory 
of elliptic curves, and can be expressed either in terms of the denominators of the 
co-ordinates of the multiples of a particular point on a particular elliptic curve, or 
in terms of special values of certain Jacobi theta functions.' 

Below, I detail the continued fraction expansion of the square root of a monic 
polynomials of even degree. In the quartic and sextic cases I observe explicitly 
that the parameters appearing in the expansion yield integer sequences defined by 
relations including and generalising that of the example Q). 

However it is well known, see for example Adams and Razar Q, that each step 
of the continued fraction expansion corresponds to addition of the divisor at in- 
finity on the relevant elliptic or hyperelliptic curve; that readily explains Zagier's 
explanation. 

1. Some Brief Reminders 
1.1. The numerical case 

We need little more than the following. Suppose w is a quadratic irrational integer 
oj , defined by lo 2 — tui + n = , and greater than the other root ZU of its defining 
equation. Denote by a the integer part of u> and set Pq = a — t, Qo = 1. Then 
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the continued fraction expansion of loq := (lo + Pq)/Qo is a two-sided sequence of 
lines, h in Z, 

lo + P h TD + P h +i . . 

—p^ = a h ; m brief uj h = a h - p h , 

Wh Wh 

with [u + Ph + i)(p + Ph+i) = —QhQh+i defining the integer sequences (Ph) and 
(Qh)- Obviously Qo divides (lo + Pq)(uJ + Pq) . This suffices to ensure that the 
integrality of the sequence (a^) of partial quotients guarantees that always Qh 
divides the norm n + tPh + P% = (^ + Ph)(w + Ph) ■ 

Comment 1. Consider the Z-module i^ ~ (Qh, Ph + lo) . It is a pleasant exercise 
to confirm that ih is an ideal of the domain Z[w] if and only if Qh does indeed 
divide the norm (lo + Ph)(Lo + Ph) of its numerator. 

One says that a real quadratic irrational LOh is reduced if 

u>h > 1 and its conjugate cJh asatisfies — 1 < i~Oh < 0. 

If the partial quotient ah is always chosen as the integer part of u>h then loq reduced 
entails all the u>h and ph are reduced; and, remarkably, ah — which starts life as 
the integer part of ujh — always also is the integer part of ph ■ Then conjugation 
of the continued fraction tableau retrieves the negative half of the expansion of loq 
from the continued fraction expansion of po . 

Comment 2. What a continued fraction expansion does. Suppose a = [ao , a% , , . . . ] 
with a > 1 , is our very favourite expansion, so much so that sometimes we go quite 
alpha — expanding arbitrary complex numbers (3 = (3q by the 'alpha' rule 

Ph = a h + (Ph - a h ) and ((3 h - a-h)^ 1 = Ph+l = o-h+i H etc. 

What can one say about those 'alpha d' complete quotients f3h ? 

Quite a while ago, in 1836, Vincent reports that either (i) all /3ft, > 1, in which 
case (3 = a; or (ii) for all sufficiently large h, \(3h\ < 1 and the real part 3?/3^ of 
Ph satisfies — 1 < diPh < 0; in other words, all those Ph lie in the left hand half of 
the unit circle. Proof: Straightforward exercise; or see [4j or, better, the survey [2j. 

Of this result Uspensky |2J] writes: 'This remarkable theorem was published by 
Vincent in 1836 in the first issue of Liouville's Journal, but later [was] so completely 
forgotten that no mention of it is found even in such a capital work as the Enzy- 
clopadie der mathematischen Wissenschaften. Yet Vincent's theorem is the basis of 
the very efficient method for separating real roots ... '. 

The bottom line is this: When we expand a real quadratic irrational a then, 
willy-nilly, by conjugation we also expand a. By Vincent's theorem, its complete 
quotients eventually arrive in the left hand- half of the unit circle and, once 'reduced', 
they stay that way. 

One readily confirms that the integers Ph and Qh are bounded by 

< 2P h + t < to - ZJ and < Qh < lo - lo . 

It follows by the box principle that the continued fraction expansion of lo is periodic. 
More, the adjustment whereby we replace lo by loq = lo + a — t arranges that loq 
is reduced. Yet more, by conjugating the tableau one sees immediately that (an 
observation credited to Galois) for any h the expansion of LOh is purely periodic. 

Comment 3. On conjugating the tableau, a putative preperiod becomes a 'post- 
period' in the expansion of p^ ; which is absurd. 
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1.2. The function field case 

Here, 'polynomial' replaces 'integer'. Specifically, set Y 2 = D(X) where D is a 
monic polynomial of even degree cleg D = 2g + 2 and defined over the base field, 
K say. Then we may write 

D(X) = A(X) 2 + AR(X) (2) 

where A is the polynomial part of the square root Y of D , so deg A = g + 1 , 
and the remainder R satisfies deg R < g. It is in fact appropriate to study the 
expansion of Z := ~(Y + A). Plainly 

C : Z 2 - AZ - R = with deg Z = g + 1 and deg Z < . (3) 

Comment 4. Note that Y is given by a Laurent series A + d-iX^ 1 + d^X^ 2 + 
■■■ , an element of K((X -1 )) ; in effect we expand around infinity. Felicitously, 
by restricting our attention to Z , and forgetting our opening remarks, the story 
we tell below makes sense over all base fields of arbitrary characteristic, including 
characteristic two. However, for convenience, below we mostly speak as if IK = Q. 

In the present context we study the continued fraction expansion of an element 
Zq of K(Z, X) leading to the expansion consisting of a tableau of lines, h £ Z, 

Z + Ph Z + P h+1 . . - 

— 7t = afl 7t ' m briet Zh = ah ~ ^ h ' sa y' W 

Wh Wh+l 

initiated by the conditions degPo < g, degQo < g and Qo divides the norm 
(Z + Pq)(Z + Pq) = —R + Pq(A + Pq) . Indeed, the story is mutatis mautandis 
precisely as in the numerical case, up to the fact that a function of K(Z, X) is 
reduced exactly when it has positive degree but its conjugate has negative degree. 
Here, analogously, we find that therefore all the and Qh satisfy 

deg Ph < g and deg Q h < g , (5) 

the conditions equivalent to the Z/> and R^ all being reduced. 



1.2.1. Quasi-periodicity 

If the base field K is infinite then the box principle does not entail periodicity. In 
a detailed reminder exposition on continued fractions in quadratic function fields 
at Section 4 of we are reminded that periodicity entails the existence of a 
non-trivial unit, of degree m say, in K[Z, X] . Conversely however, the exceptional 
existence of such a unit implies only 'quasi-periodicity' — in effect, periodicity 
'twisted' by multiplication of the period by a nonzero element of K. The existence 
of an exceptional unit entails the divisor at infinity on the curve C being torsion 
of order dividing m. If quasi-periodic, the expansion of the reduced element Zq = 
(Z + Pq)/Qo is purely quasi-periodic. 

Comment 5. Consider the surprising integral 

/6a: dx 
Vx 4 + Ax 3 - 6x 2 + 4.x + 1 
= log(a; 6 + 12a; 5 + 45a; 4 + 44a; 3 - 33a; 2 + 43 

+ (x 4 + 10a; 3 + 30x 2 + Tlx - 11) x 4 + 4jj 3 - 6a; 2 + 4a; + 1 ) , 
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a nice example of a class of pseudo-elliptic integrals 



/ 



f{x)dx 



log(a{x) + b(x)y/D(x)). 



(6) 



Here we take D to be a monic polynomial denned over Q , of even degree 2g + 2 , 
and not the square of a polynomial; /, a, and b denote appropriate polynomials. 
We suppose a to be nonzero, say of degree m at least g + 1 . One sees readily that 
necessarily deg b — m — g — 1 , that deg / = g , and that / has leading coefficient 
to. In the example, to = 6 and g = 1. 

The trick is to recognise that obviously a 2 — b 2 D is a non-zero constant and 
a + by/15 is a unit of degree to in the domain Q(x, yj D(x) and is not necessarily 
of norm ±1 — it is this that corresponds to gwosj-periodicity; for details see |l4|. 

1.2.2. Normaility of the expansion 

In the sequel, I suppose that Zq has been so chosen that its continued fraction 
expansion is normal: namely, all its partial quotients are of degree 1. This is the 
generic case if IK is infinite. Since I have the case K = Q in mind, I refer to elements 
of K as 'rational numbers'. 

Comment 6. Naively, it is not quite obvious that the case all partial quotients of 
degree one is generic, let alone that this generic situation can be freely arranged 
in our partcular situation. I comment on the latter matter immediately below but 
the former point is this: when one inverts d-\X + d-2X + • • ■ , obtaining a 
polynomial plus e„iX _1 + e_2^"~ 2 + ■ • ■ it is highly improbable that e_i vanishes 
because e_i is a somewhat complicated rational function of several of the d-i . See 
my remarks in |l5| . 

1.2.3. Ideal classes 

If K is infinite, choosing Pq and Qo is a matter of selecting one of infinitely many 
different ideal classes of Z-modulcs {(Qh> Z + Ph) '■ h € Z}. Only a thin subset of 
such classes fails to give rise to a normal continued fraction expansion. Of course 
our choice of Pq and Qo will certainly avoid the blatantly singular principal class: 
containing the ideal (1, Z) . 

1.2. 4- Addition on the Jacobian 

Here's what the continued fraction expansion does. The set of zeros {uJh.i, ■ ■ ■ , ^h,g} 
of Qh defines a rational divisor on the hyperelliptic curve C . In plain language: any 
one of these zeros u>h is the X co-ordinate of a point on C , here viewing C as 
defined over the algebraic extension K^^) , generically of degree g over the base 
field IK. In particular, in the case g — 1 , when C is an elliptic curve, the unique zero 
Wh of Qh provides a point on C defined over the base field; see H2.2I at page 12171 
below. 

Selecting Pq and Qo is to choose a divisor class, say M, thus a point on the 
Jacobian Jac(C) of the curve C defining Z. Let S denote the divisor class defined 
by the divisor at infinity. Then each complete quotient Zh = [Z + Ph)/Qh has 
divisor in the class Mh+i '■= M + hS . I show this in 01 forg — 1 making explicit 
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remarks of Adams and Razar See further comment at H2.2I For higher genus 
cases, one wilkl find helpful the introduction to David Cantor's paper |j| and the 
instructive discussion by Kristin Lauter in . A central theme of the paper Q is 
a generalisation of the phenomenon to Pade approximation in arbitrary algebraic 
function fields. My suggestion that partial quotients of degree one are generic in 
our examples is the same remark as that divisors on C/K typically are given by g 
points on C/F , where F is some algebraic extension of K . 



2. The Continued Fraction Expansion 

Evidently, the polynomials Ph and Qh in 

z + Ph =ah _ Z + P h+1 

Qh Qh+1 

are given sequentially by the formulas 

Ph + Ph+i + A = a h Q h and 

-QhQh+i = {Z + Ph+i)(Z + Ph+i) = -R + P h +i{A + P h+1 ). (8) 



2.1. A naive approach 



At first glance one might well be tempted to use this data by spelling out the 
first recursion in terms of g equations linearly relating the coeficients of the Ph , 
and the second recursion in terms of 2g + 1 equations quadratically relating the 
coefficients of the Qh- Even for g = 1, doing that leads to a fairly complicated 
analysis; see 0] , and the rather less clumsy 01 • F°r .9 = 2,1 had to suffice myself 
with the special case degi? = 1 so that several miracles (I called it 'a ridiculous 
computation') could yield a satisfying result I realised that a quite different 
view of the problem would be needed to say anything useful in higher genus cases. 



2. 2. Less explicit use of the recursions 

Denote by cuh a zero of Qh{X) . Then 

Ph(uj h ) + Ph+iM + A(w h ) = 

so - QhQh+i = -R+ Ph+i (A + P h+ i) 

becomes R(uj h ) = -Ph+i(oJh)Ph(uh)- (9) 

This, together with a cute trick, already suffices to tame the g = 1 case: that is, 
the case beginning as the square root of a monic quartic polynomial. Indeed, if 
g = 1 then degP^ = g — 1 ; so the Ph are constants, say -P/i(A) —: eh- Also, 
degQ/j = 1, say Q h {X) =: v h (X - w h ), so ui h = w h . Further, degR < 1, say 
R(X) =: v(X-w). 

First, © tells us directly that v(w — Wh) = ehCh+i- Second, this is the 'cute 
trick', obviously —R(wh) = v(w — Wh) = Qh(w) ■ v/vh- An intelligent glance at |(HJ 
reminds us that Vh-iVh — —eh - Hence, by JSJ and R(w) = 0, 



e h ^xe\e h+ i = -Qh-i{w)Q h {w) ■ v 2 /e h = v 2 (e h + A(w)) 
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The exciting thing about the recursion 

e h -ie\e h +i = v 2 (e h + A(w)) (10) 

is that, among the parameters varying with h, it involves the eh alone; moreover, 
its coefficients v 2 and v 2 A(w) depend only on the curve C and not on the initial 
conditions Pq = e$ and Qq = vq(X — wq) ■ 

More, if X = Wh then Z = —e% or Z = —eh+i; yielding pairs of rational points 
on C . Moreover, the transformation U = Z , V — v = XZ transforms the curve C 
to a familiar cubic model 

£ : V 2 — vV — monic cubic in U with zero constant coefficient, (11) 

essentially by moving one of the two points Sc , say, at infinity on C to the origin 
Ss = (0,0) on £. As at i jl.Lj.4l on page 12 1^1 above, denote the point (wh,—eh) on 
C by Mh+\ ■ Recall that an elliptic curve is an abelian group with group operation 
denoted by +. Set M\ — M. One confirms, see 0] for details, that Mh+i = 
M + hSc by seeing that this plainly holds on £ where the addition law is that 
F + G + H = if the three points F , G , H on £ lie on a straight line. 

I assert at m. 2. 41 that precisely this property holds also in higher genus. There, 
however, one is forced to use the gobbledegook language of 'divisor classes on the 
Jacobian of the curve' in place of the innocent 'points on the curve' allowed in the 
elliptic case. 

Comment 7. By the way, because each rational number —eh is the U co-ordinate 
of a rational point on £ it follows that its denominator must be the square of an 
integer. 

2. 3. More surely useful formulas 

Set a,h(X) = (X + Vh)/uh \ so Uh is the leading coefficient of Qh (its coefficient 
of X 9 ) and ah vanishes at —Vh- Below, we presume that dh denotes the leading 
coefficient of Ph (its coefficient of X 9-1 ). Then 

R(X) + Ph+i(X)Ph(X) = Qh(X)Q h +i(X) + P h+1 (X)(X + is h )Qh(X)/u h 

= Qh-i(X)Q h (X) + P h (X){X + Vh )Qh{X)/u h . 

Plainly, we should divide by Qh(X)/uh and may set 

C h {X) := (R{X) + Ph +1 (X)Ph(X))/(Qh{X)/uh) 
= u h Q h+ i{X) + P h+1 {X)(X + v h ) 

= UhQh-i(X) + Ph{X){X + is h ). (12) 

Since deg R < g , and the P all have degree g — 1 , and the Q degree g , it follows 
that the polynomial C has degree g — 2 if g > 2 , or is constant in the case g = 1 . 
If g = 2 it of course also is constant and then, with R(X) = u(X 2 — vX + w) , its 
leading coefficient is dhdh+i + u so we have, identically, Ch{X) = dhdh+i + u. If 
g > 3 then Ch is a polynomial with leading coefficient dhdh+i ■ 

2.4. The case g — 2 

If deg A = 3 then R = u(X 2 — vX + w) is the most general remainder. In the 
continued fraction expansion we may take Ph(X) — dh(X + eh)- As just above, 
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denote by Uh the leading coefficient of Qh{X) and note that Uh-iUh = —dh- In 
the case g — 1 we succeeded in finding an identity in just the parameters eh ; here 
we seek an identity just in the dh ■ 

First we note, using ljl2|l and Ch(X) = dhdh+i + u, that 

dhdh+i + u = C h (-e h ) = u h R(-e h )/Qh(-e h ) = u h Q h -x(-e h ) ■ (13) 
dh-A +u = C h -i(-e h ) = u h -iR{-e h )/Qh-x{-eh) = u h -iQh{-eh) ■ (14) 

Hence, cutely, 

(dh-idh + u)(d h dh+i + u) = Uh-iUhR(-e h ) = -ud h {e\ + ve h + w) . (15) 
Set R(X) =: u(X ~ uo)(X -W). Also from O, we nav e 

C h (uj)Q h (uj) = P h +i{u)Ph{u) 

and therefore 

Ch~i(u)C h (to)Qh-i(u)Qh(uj) = P h -!{u;)(P h (u)) 2 P h+1 (u;). 
But -Q h -i(u)Qh(u) = Ph(u){A(u) + P h {u)) ■ Together with CE3) we find that 
u 3 (A(uj) +d h (oj + e h ))(A(uJ) + d h (U + e h )) 

= d h -xd\dh+\{dh-2dh-x + u)(d h+ idh+2 + u). (16) 

In the special case u = 0, that is: R(X) = —v(X — w),& little less argument yields 
the more amenable 

-v 3 (A(uj) +d h (w + e h )) = d h ^ 2 dl-idldl +1 d h+ 2 ■ 

In this case, we have —vdh(w + eh) = —dhR(-eh) = dh-id^dh+i , finally providing 

d h -2dl_idhdl +1 dh+2 = v 2 d h -\d\d h+1 - v 3 A(w) . (17) 

Comment 8. This reasonably straightforward argument removes the 'miraculous' 
aspects from the corresponding discussion in [l^. Moreover, it gives a result for 
ii ^ 0. However, I do not yet see how to remove the dependence on in 11611 so 
as to obtain a polynomial relation in the ds. 

3. Somos Sequences 

The complexity of the various parameters in the continued fraction expansions 
increases at frantic pace with h . For instance the logarithm of the denominators 
of eh of the elliptic case at ij2.2l is readily proved to be 0(h 2 ) and the same must 
therefore hold for the logarithmic height of each of the parameters. 

Denote by A\ the denominator of eh (recall the remark at Comment [7] on 
page 121 71 that these denominators in fact are squares of integers). It is a remarkable 
fact holding for the co-ordinates of multiples of a point on an elliptic curve that in 
general 

Ah-iA ll+ i = e h A\ (18) 

Moreover, it is a simple exercise to see that Ijl8|l entails Ah- 2 Ah+2 = &h~\&h£h+\A\ . 
Thus, on multiplying the identity I|1Q(I . namely eh-\e\eh+\ = v 2 (e,h + A{wy\ , by 
A\ we find that 

Ah- 2 A h+ 2 = v 2 A h -iA h+1 + v 2 A(w)A 2 (19) 
gives a quadratic recursion for the 'denominators' Ah ■ 
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Comment 9. My remark concerning the co-ordinates of multiples of a point on an 
elliptic curve is made explicit in Rachel Shipsey's thesis |2l| . The fact this also holds 
after a translation is shown by Christine Swart [2^j | and is in effect again proved 
here by way of the recursion l|10|l . My weaselling 'in general' is to avoid my having 
to chat on about exceptional primes — made evident by the equation defining the 
elliptic curve — at which the Ah may not be integral at all. In other words, in 
true generality, there may be a finite set S of primes so that the Ah actually 
are just 5 -integers: that is, they may have denominators but primes dividing such 
denominators must belong to S . 

3.1. Michael Somos' 's sequences 

Some fifteen years ago, Michael Somos noticed 0,H3], that the two-sided sequence 
Ch-2Ch+2 = G'h-iCh+i+Cfr , which I refer to as 4-Somos in his honour, apparently 
takes only integer values if we start from Ch-x = Ch = Ch+i = Ch+2 = 1. 
Indeed Somos went on to investigate also the width 5 sequence, Bh-2Bh+3 — 
Bh-iBh+2 + Bh.Bh+i , now with five initial Is, the width 6 sequence Dh-sDh+z = 
Dh—iDh+i l-Dfo+i +Df t , and so on, testing whether each when initiated by an 

appropriate number of Is yields only integers. Naturally, he asks: "What is going 
on here?" 

By the way, while 4-Somos (A006720) 5-Somos (A006721), 6-Somos (A006722), 
and 7-Somos (A006723) all do yield only integers; 8 -Somos does not. The codes in 
parentheses refer to Neil Sloane's 'On-line encyclopedia of integer sequences' [22$ . 

3.2. Elliptic divisibility sequences 

Sequences generalising those considered by Somos were known in the literature. 
Morgan Ward had studied anti-symmetric sequences (Wh) satisfying relations 

W h . m W h+m Wl = Wh- n Wh +n Wl - W m . n W m+n Wl . (20) 

He shows that if W\ — 1 and \V2\W4 then 6 implies that W W&; that is, the 
sequences become divisibility sequences (compare the Fibonacci numbers). For a 
brief introduction see Chapter 12 of 0. 

There is a drama here. The recurrence relation 

= WlWh-tWh+t - w x w 3 wl , 

and four nonzero initial values, already suffices to produce (Wh) ■ Thus H2()|l for all 
m and n is apparently entailed by its special case n = 1 and m = 2 . The issue is 
whether the definition l|2U|) is coherent. 

One has to go deep into Ward's memoir j2j| to find an uncompelling proof that 
there is in fact a solution sequence, namely one defined in terms of quotients of 
Weierstrass sigma functions. More to the point, given integers W\ = 1 , W2 , W3 , 
and W4 , there always is an associated elliptic curve. In our terms, there is a curve 
C : Z 2 — AZ — R = with deg^4 = 2, degi? = 1 and the sequence (Wh) arises 
from the continued fraction expansion of Z\ = Z/(—R). I call (Wh) the singular 
sequence because in that case e\ = — so the partial quotient (Iq(X) is not linear. 

My 'translated' sequences (Ah) were extensively studied by Christine Swart in 
her thesis E^ ]. 
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3.3. Somos sequences 

It is natural to generalise Michael Somos 's questions and to study recurrences of the 
kind he considers but with coefficients and initial values not all necessarily 1 . Then 
our recurrence Q19f) is a general instance of a Somos 4 sequence, an easy computation 
confirms that — v = W2 and v 2 A(w) = — W3, and given the recursion and four 
consecutive Ah one can readily identify the curve C and the initial 'translation' 
M = (w , -e ). 

Comment 10. One might worry (I did worry) that Ah-2Ah+2 = aAh-iAh+i+bA^ 
does not give a rationally defined elliptic curve if a is not a square. No worries. 
One perfectly happily gets a quadratic twist by a of a rationally defined curve. 

For example, 4-Somos, the sequence (Ch) = (..., 2, 1, 1, 1, 1, 2, 3, 7, ... ) 
with C h - 2 C h+2 = C h -iC h+ i + C\ arises from 

C : Z 2 - (X 2 - 3)Z - (X - 2) = with M = (1, -1) ; 

equivalently from £ : V 2 - V = U 3 + 3U 2 + 2U with M £ = (-1, 1) . 

Christine Swart and I found a nice inductive proof 0] that if (Ah) satisfies (|19|l 
then for all integers m and n, 

A h - m A h+m W 2 = W 2 n A h - n Ah+ n - W m - n W m+n A 2 . 

Our argument obviates any need for talk of transcendental functions and is purely 
algebraic. 

It also is plain that A h -iA h+ i = e h Al yields A h -iA h+ 2 = e h eh+iA h Ah+i and 
Ah-^Ah+s = e h ~ie 2 l e 2 h+1 eh+2AhA h+ \ . However, although JTUJl directly entails 

e h -ie 2 h e 2 h+l e h+ 2 = -v 2 A(w)e h e h +i + (v 4 + 2wv 3 A(w)), 

it requires some effort to sec this. Whatever, (|19ll eventually also gives 

WiW 2 A h ^ m A h+m+1 = W m W m+1 A h -iA h+2 - W m -xW m +2A h A h+ i . (21) 

For details see [13] and |lflj . 

The case m — 2 of 1121(1 includes all Somos 5 sequences. The sequence 5 -Somos, 
(B h ) = (..., 2, 1, 1, 1, 1, 1, 2, 3, 5, 11, ... ) with Bh- 2 B h+3 = B h -iB h+2 + B h B h+1 , 
arises from 

Z 2 - (X 2 - 29)Z + 48(X + 5) = with M = (-3, -8) ; 

equivalently from £ : V 2 + UV + 6V = U 3 + 7U 2 + 12U with M £ = (-2, -2) . 

Actually, a Somos 5 sequence (Ah) may also be viewed as a pair (A2h) and 
(A 2 h+i) of Somos 4 sequences coming from the same elliptic curve but with different 
translations (in fact differing by half its point S ) ; see the discussion in • 

3.4- Higher genus Somos sequences 

My purpose in studying the elliptic case was to be able to make impact on higher 
genus cases. That's been only partly achieved, what with little more than 117fl at 
page !218l to show for the effort. I tame (|T7jl by defining a sequence (Th) , one hopes 
of integers, by way of Th-iTh+i = dhT 2 . We already know that then Th—^Th+n = 
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dh-xd^dh+xT^; also T h ^ 3 T h+3 = d h - 2 d h \_ l d\d 2 l+l d h +2Tl follows with only small 
extra effort. Thus (|17|l yields 

T h - 3 T h+3 = v 3 T h _ 2 T h+2 - v 2 A(w)T 2 . (22) 

Then the sequence (T h ) = (..., 2, 1, 1, 1, 1, 1, 1, 2, 3, 4, 8, 17, 50, . . .) 
satisfying Th- 3 T} l+3 = T^^Th+i + Tj 2 is readily seen to derive from the genus 2 
curve 

C : Z 2 - Z(X 3 -AX + 1) - (X - 2) = 0. (23) 
A relevant piece of the associated continued fraction expansion is 

Z + 2X -1 Z + X 



X 2 -l X 2 -l 

Z+X Z + X-l 

-JC 



-(X 2 -2) -{X 2 ~2) 
Zo:=j +X -\=X + l * + X - 



X 2 - X - 1 X 2 - X 

Z+X-l Z+X 
A — 



-(X 2 -2) -{X 2 -2) 

Z + X Z + 2X-1 

X — 



X 2 - 1 X 2 - 1 



illustrating that M is the divisor class defined by the pair of points {<p,~0) and 
(jp, ip) — here, <p is the golden ratio, a happenstance that I expect will please 
adherents to the cult of Fibonacci. 

Comment 11. There does remain an issue here. Although i|l(J|) is concocted on 
the presumption that is never zero, it continues to make sense if some 
should vanish — for instance in the singular case. However, l|17|) is flat out false if 
dh-xdhdh+i = 0. Indeed, a calm study of the argument yielding l(T7|l sees us divid- 
ing by v(eh + w) at a critical point. These considerations together with Cantor's 
results |8| suggest that should always be reported as multiplied by Th and, 
more to the point, that the general genus 2 relation when ii / will be cubic 
rather than quadratic. 



4. Other Viewpoints 

My emphasis here has been on continued fraction expansions producing sequences 
M + hS of divisors — in effect the polynomials Qh — obtained by repeatedly 
adding a divisor S to a starting 'translation' M . That viewpoint hints at arithmetic 
reasons for the integrality of the Somos sequences but does not do that altogether 
convincingly in genus greater than one. 



4-1. The Laurent phenomenon 

As it happens, the integrality of the Somos sequences is largely a combinatorial phe- 
nomenon. In brief, as an application of their theory of cluster algebras, Fomin and 
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Zelevinsky |l(j prove results amply including the following. Suppose the sequence 
(yh) is defined by a recursion 

Uh+nUh = ayh+rljh+n-r + ^Vh+sUh+n-s + JUh+tUh+n-t , 

with < t < s < r < ^n. Then the yh are Laurent polynomials* in the variables 
yo , 2/1 , ■ ■ • , y n -i and with coefficients in the ring Z[a, (3, 7] . That deals with all 
four term and three term quadratic recursions and thus with the cases Somos 4 to 
Somos 7. Rather more is true than may be suggested by the given example. 

4-2. Dynamic methods 

Suppose we start from a Somos 4 relation Ah~2^h+2 = ctAh-iAh+i + flA^ 
and appropriate initial values A , Ai, A 2 , A 3 . Then one obtains rationals eft = 
Ah-i Ah+i/A\ satisfying the difference equation 

e-h+i = [d H J . 

eh-ieh v eft/ 

The point is that this equation has a first integral given by 

J := J(e h -i, eft) = eft_ieft + a( 1 ) H — = J(e h , e h+ i) 

\eh-i eft/ eft_ieft 

and one can now construct an underlying Weierstrass elliptic function p. Indeed, 
the readily checked assertion that given y £ C there are constants a and so that 

(p(x + y) - p(y))(p(x) - p{yj) (p(x - y) - p{y)) = -a(p(x) - p(y)) +(3 

reveals all; particularly that a — p'(y) 2 , (3 — p'(y) 2 (p(2y) — p(y)) . 
Specifically, after fixing x by J = p" (x) one sees that 

-eft = p[x + ny) - p(y) . 

A program of this genre is elegantly carried out by Andy Hone in |l2| for Somos 4 
sequences. In 0, the ideas of 01 are shown to allow a Somos 8 recursion to 
be associated with adding a divisor on a genus 2 curve. Incidentally, that result 
coheres with the guess mooted at Comment ^2 above that there always then is a 
cubic relation of width 6 . 
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